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We study families of solitons in a two-dimensional (2D) model of the light transmission through 



o . 

, a photorefractive medium equipped with a (quasi-)one-dimensional photonic lattice. The soliton 



families are bounded from below by finite minimum values of the peak and total power. Narrow 
solitons have a single maximum, while broader ones feature side lobes. Stability of the solitons is 
checked by direct simulations. The solitons can be set in motion across the lattice (actually, made 
' tilted in the spatial domain), provided that the respective boost parameter does not exceed a critical 

value. Collisions between moving solitons are studied too. Collisions destroy the solitons, unless 
^ their velocities are sufficiently small. In the latter case, the colliding solitons merge into a single 

' stable pulse. 

d ' 

PACS numbers: 



I. INTRODUCTION 
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' As was first predicted in Refs. a periodic lattice potential induced in a photorefractive medium, which is 
' characterized by the saturable nonlinearity, may be an efficient tool for creation and stabilization of two-dimensional 



[ spatial solitons of various types. The prediction was followed by the experimental observation of ordinary solitons 

locahzed vortices with the topological charge 1 belonging to the first (lowest) Q or second Q 



o 

. [21, localized vortices with the topological charge 1 belonging to the first (lowest) [3j or second bandgap in the 
respective linear spectrum (stable higher-order vortices and supervortices in such systems were recently predicted too 
' Q), dipole- and quadrupole-mode solitons Q|, steady patterns in the form of soliton necklaces 7], and some others. 
^ ■ A review of the field was recently given in Ref . Q . 

' In the experiment, the periodic potential is induced by the photonic lattice, which is created as a superposition 
CZ of counterpropagating laser beams illuminating the photorefractive crystal in the ordinary polarization, in which the 
light propagation is nearly linear; then, the spatial solitons are built in a probe beam launched through the lattice 
in the extraordinary polarization (i.e., polarized along the crystalline c axis), which is subject to strong nonlinearity 
induced by the dc bias electric field applied to the crystal QQ. Using this technique, the square photonic lattice can 
be created by two pairs of counterpropagating beams illuminating the crystal in directions orthogonal to each other 
and to the probe beam that gives rise to the soliton(s). 

On the other hand, 2D solitons can also be supported by a low-dimensional, i.e., quasi-one-dimensional (QID) 
periodic potential, that may be readily induced by a single pair of counterpropagating beams illuminating the bulk 
crystal in the ordinary polarization. Recently, it has been predicted that the QID lattice may efficiently stabilize 2D 
solitons in the model with the cubic (rather than saturable) nonlinearity jsj. This result directly applies to Bose- 
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Einstein condensates with attractive interactions between atoms, trapped in a photonic lattice; moreover, it has also 
been demonstrated that a quasi-2D lattice can stabilize 3D solitons against strong collapse that the cubic 

nonlinearity gives rise to in the latter case. 2D solitons supported by the QID lattice naturally demonstrate strong 
anisotropy, which makes them essentially different from the usual 2D solitons. A significant advantage offered by the 
use of the QID lattice is the fact that the remaining free direction allows the solitons to move (in the spatial domain, 
"motion" means a tilt of the soliton beam) , which opens a way to study collisions between them, formation of bound 
states, etc. 9']. The mobility of 2D solitons may also be strongly anisotropic in some 2D lattices 

The objective of the present work is to introduce 2D solitons in the model of the photorefractive media with the 
saturable nonlinearity and QID periodic potential (QID solitons - which are, effectively, one-dimensional objects - 
in the model of the 2D photorefractive medium with the QID lattice were introduced in Ref. J^)- The soliton 
solutions are constructed in Section 2. We demonstrate that the QID solitons are characterized by minimum peak 
and total intensities necessary for their existence, and they are stable in the entire existence region. Moving solitons 
and collisions between them are studied in Section 3. We find that collisions destroy the solitons, unless the collision 
"velocity" (in fact, the relative tilt of the two spatial solitons) is small enough; in the latter case, the colliding solitons 
merge into a single one, irrespective of the orientation of the velocity vector relative to the QID lattice. 

The model outlined above is based on an equation for the spatial evolution of the probe field (the slowly varying 
amplitude U of the electromagnetic wave in the extraordinary polarization) , which follows the standard description 
of photorefractive media Jj. In normalized units, the equation takes the form: 



dU d^U d^U U ^ 

l--^ ^ -TT^ + -TT^ 9 = 0, (1) 

dz dx^ dy^ l + /oCos2(7ra;/d) + |C/|' 

where /q and d are the peak intensity and period of the photonic lattice induced by the superposition of counterprop- 
agating waves launched along the x axis in the ordinary polarization. The normalized propagation and transverse 
coordinates, z and {x,y), are proportional to their counterparts measured in physical units, Z and (X,Y), so that 
z = Z / (koAno) and {x,y) = {X,Y) / ^/2konoAno, where fco is the propagation constant of the probe wave, and 
Ano = nQr33E/2 is the change of the refractive index no (accounted for by the electro-optic coefhcient raa) caused by 
the dc electric field E. 



II. TWO-DIMENSIONAL SOLITONS AND THEIR STABILITY 



Soliton solutions to Eq. ^ are searched for as 



C/ = M(a;,y)e-'^^ 

where — /i is the shift of the propagation constant in the soliton, and the real function u satisfies the equation 



dx^ dy^ 



1 



u = 0. 



(2) 



(3) 



1 + /o cos^ {nx/d) + V? 

Solutions u{x, y) of this equation were found by means of an iteration procedure in the Fourier space. To this end, 
following a numerical method elaborated in Ref. 6], Eq. (j2Jl was rewritten for the Fourier transform u{kj:, ky) in the 
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form 

u=—-^ J^f /oCOS^(Wrf)-. \vlu \\ 

kl + kl-^i+\\ \^l + /ocos2(7rx/d) + |u|7 \^l + /ocos2(^x/d) + |u|7 J 

standing for the 2D Fourier transform. A direct iteration procedure applied to Eq. ^ does not converge, in 
the general case. Therefore, the equation was modified as follows: defining integral factors, 

k^, — U + ^^u — T \ ^ ttU I > u*dkj:dky (5) 

" ' \l + /ocos2(7rx/(i) + |u|' J j 

(3 ^ I IT\ — J u*dk^dky^ (6) 

/o cos"^ [nx/dj + \u\ ) 

where * stands for the complex conjugation, the following iterative equation was introduced. 





/q COS^ [TTx/d) + \Un\^ J 




(7) 



/q cos^ {ttx /d) 

where a„ and are the factors (01 and jnj corresponding to the function m„. Fixed points of Eq. corresponding 
to lim (a„//3„) — 1, yield solutions to Eq. as well. The iterative procedure based on Eq. (0) provides for fast 

n ^oo 

convergence, and produces solutions displayed below. 

Typical examples of the 2D solitons are shown in Fig. ^ In particular, the picture observed in the left column of 
the figure is typical to cases when the lattice potential is weak, and/or the soliton's peak intensity essentially exceeds 
the lattice's amplitude /q: the soliton is practically isotropic, without a conspicuous effect of the lattice. The picture 
in the right column is typical for a relatively strong lattice with a large period: then, the soliton is almost entirely 
trapped in one potential trough, assuming an elliptic shape. In either case, the soliton's shape features a single lobe, 
either circular or elliptic one. 

The most interesting typical case is represented by the central column in Fig. ^ which demonstrates a soliton with 
well-pronounced side lobes in a moderately strong lattice (in physical units, this case corresponds to typical values 
of parameters available in the experiment). Naturally, the multi- humped structure is observed only along the axis x, 
while in the free direction, y, the soliton always features a simple single-hump form. Results reported below are given 
for the same period, d = 2tt, which gives rise to the central column in Fig. ^ For other values of d which are neither 
very small nor very large, the results are quite similar. 

The solitons are characterized by the peak intensity (power), Ip = \u[x — 0,y — 0)|^, and integral intensity, 

P + CO P + CO 

P= / \u{x,y)fdxdy. (8) 



Accordingly, soliton families for given values of Iq and d are represented by the dependences P(/i) and Ip{fJ,), see Fig. 
121 [recall — /i is the soliton propagation-constant shift, defined in Eq. ©J. An important feature observed in this figure 
is that the 2D solitons exist, with a given lattice strength /q, only for P and Ip exceeding certain finite minimum 
(threshold) values, Pmin and /p^min; for instance, /p.min = 0.0411 for Iq — 5 and d ~ 2tt. Both minimum values are 
shown, as functions of the strength Iq of the photonic lattice, in Fig. 13 It should be noticed that a lower intensity 
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threshold necessary for the existence of 2D sohtons in models combining lattice potentials with various nonlinearities 
was found in earlier works 

It is easy to check that the soliton families shown in Figs. ^ and |21 belong to the semi-infinite bandgap in the 
spectrum of the linearized equation It is known that, on top of the semi-infinite gap, a sufficiently strong lattice 
may give rise to solitons in finite bandgaps of models combining a periodic potential and saturable nonlinearity (see, 
e.g., Ref. jUij). In this work, we do not aim to look for gap solitons of such types. 

Proceeding to the investigation of stability of the 2D solitons, we first of all note that the negative slope of the de- 
pendence ^'(/J.), obvious in Fig. 12 suggests possible stability of the soliton families pursuant to the Vakhitov-Kolokolov 
(VK) criterion 15], which may guarantee the absence of unstable modes of small perturbations corresponding to real 
eigenvalues of the instability growth rate. The full stability is not provided by the VK criterion, and, moreover, the 
applicability of the criterion even to perturbations with real eigenvalues was not proven in the present context. In fact, 
a counter-example is known: soliton subfamilies which should be VK-unstable in a ID model combining a periodic 
potential of the Kronig-Penney type and cubic-quintic nonlinearity (another variety of the saturation) were found to 
be completely stable |l4| . 

We have tested stability of the 2D solitons in direct simulations of Eq. , making use of the fast Fourier transform 
in directions x and y and Runge-Kutta method to advance in z. In each simulation, some amount of random noise was 
added to the soliton as a perturbation. It has been concluded that all the solitons are stable, as (formally) predicted 
by the VK criterion. An example illustrating the stability of a large-amplitude 2D soliton with /o = 5 and Ip = 56.76 
is displayed in Fig. 0fa). Broad solitons, with the peak intensity close to the threshold value Ip^min, are stable too, 
although their peak intensity may slowly grow with z, as shown in Fig. EJb) for the same lattice strength, Iq = 5, as in 
Fig. EJa), but Ip = 1.1. A plausible explanation to the latter effect is that the combination of the lattice potential and 
self- focusing nonlinearity (for relatively small Ip, the saturation does not suppress the self- focusing) leads to sucking 
the perturbation wave field into the spot where the soliton's maximum is located, and the build-up of the additional 
wave field around the soliton's peak may be conspicuous, against the backdrop of the relatively small value of Ip. 



A. Moving two-dimensional solitons and their collisions 

Solutions U{x,y,z) for solitons "moving" (actually, tilted) along the free direction y can be generated from the 
"quiescent" solitons reported in the previous section, U{x,y,z), by means of the Galilean transformation, 

U {x, y,z) ^ U (x, y - 2QyZ, z) exp {iQyX - iQ^z) , (9) 

where Qy is an arbitrary boost parameter, that determines the soliton's "velocity" (tilt) 2Qy. On the other hand, 
generation of solitons tilted along the x axis is a nontrivial problem. In the ID version of Eq. without the term 
d^U/dy^, a family of tilted solitons was introduced in Ref. [lj |. 

We looked for solitons "moving" along the x axis by simulating the evolution of initial states of the form 

Uo{x,y) ^U{x,y)eyi^{iQxx) , (10) 

cf. Eq. jnj, where J7(a;, y) corresponds to a zero-velocity stationary soliton solution. The boost parameter in Eq. Hl()|l . 
Qx, was gradually increased from a run to a run, until reaching a critical (maximum) value Qmax, at which the initial 
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configuration (|10|l does not generate any soliton, but rather gets destroyed into small-amplitude waves. Results of the 
simulations are summarized in Fig. [3 in the form of plots showing Qmax as a function of the soliton's peak intensity 
Ip, for different values of the lattice strength Iq. Measurement of the established value of the average "velocity" (tilt) 
of the solitons generated by the initial condition (|10|l with Qx < Qmax produces values which are quite close to ones 
{2Qx) corresponding to the Galilean transform in the free space. Finally, applying the Galilean transformation 
to the soliton already moving along the x direction, one can generate a pulse moving in any direction (in particular, 
along the diagonal, x = y, see below). 

Once the stability limits for the moving solitons are available, the next step is to consider collisions between them. 
The outcome of the collision may depend on the magnitude and direction of the velocities, and the aiming mismatch 
(its zero value corresponds to the head-on collision). 

Numerous simulations demonstrate that the collision completely destroys both solitons, unless their velocities are 
sufficiently small. An example of the destructive collision is displayed in Fig. El In particular, in the case of Iq — 5, the 
solitons must not be set in motion by the boost with \Qx \ or \Qy \ in excess of 0.05, otherwise the collision will destroy 
them. This critical value is the same for the motion in the x and y directions, up to the accuracy of the simulations. 
Moreover, if the solitons are boosted in the diagonal direction and then collide, the same limit, \Q\ = 0.05, was 



found for the absolute value of the corresponding vectorial boost, \Q\ = ■\JQx + <9yi i-^-, the critical boost is virtually 
isotropic. 

Collisions between solitons boosted by IQj with values smaller than the critical one result in their merger into a 
single pulse. An example of the merger resulting from the collision with a finite aiming mismatch is displayed in Fig. 
□ Further, the mergers caused by head-on collisions in the diagonal or vertical (y) direction are presented in Fig. |H1 
As seen from the figures, the pulse generated by the merger performs intrinsic vibrations, but remains stable. 

The above examples displayed collisions of single-lobe solitons. Collisions between their broader counterparts, which 
feature side lobes in the x-direction, are quite similar. Figure El presents an example of the head-on collision between 
sufficiently slowly moving solitons of the latter type. The collision again leads to the merger of the pulses into a single 
one, which then performs conspicuous intrinsic vibrations, but remains a stable object. Additional simulations show 
that the soliton produced by the merger of two multi-lobe ones can also easily move across the lattice, if given a push. 

Thus we conclude that collisions between the 2D solitons are always strongly inelastic (both the destruction and 
merger of the colliding solitons are inelastic outcomes), attesting to the fact that the present model is far from any 
integrable limit, where collisions between solitons would be elastic. For comparison, we note that the collisions may 
be less inelastic in the 2D model with the QID periodic potential if the nonlinearity is cubic In that model, elastic 
collisions are possible (passage of the solitons in the case of a finite aiming mismatch, or their mutual bounce after 
the head-on collision, if the phase shift between them is tt). An additional inelastic outcome of the head-on collisions 
between in-phase 2D solitons, which was observed in the model with the cubic nonlinearity, but cannot occur if the 
nonlinearity is saturable, is collapse of the single pulse formed after the merger of two solitons (i.e., formation of a 
singularity after a finite propagation distance). 
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III. CONCLUSION 

In this paper, we have proposed a model of the two-dimensional (2D) medium with the saturable nonlinearity 
and quasi-lD lattice potential, that can be realized in the spatial domain in a bulk photorefractive crystal. The 
subject of the analysis were 2D solitons (ones belonging to the semi-infinite bandgap in the linear spectrum). It was 
demonstrated that they form families bounded from below by finite minimum values of the peak and total intensities. 
Narrow pulses feature a single maximum, while broad solitons contain side lobes. Direct simulations confirm that the 
solitons are stable. They can be set in motion (actually, tilted in the spatial domain) in an obvious way along the 
quasi-lD lattice, and also across the lattice, provided that, in the latter case, the boost parameter does not exceed a 
critical value, beyond which the soliton is destroyed. 

Collisions between stable moving solitons were studied in detail, with the conclusion that the collisions destroy 
the solitons, unless their velocities are sufficiently small. In the latter case, the colliding solitons merge into a single 
stable pulse, that performs intrinsic vibrations. The predictions reported in this paper can be readily implemented 
experimentally in photorefractive crystals. 
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Fig. 1: Typical examples of two-dimensional solitons found from Eq. © with a moderately strong quasi-one-dimensional photonic 
lattice, corresponding to lo — 5. Panels (a), (b), and (c) show the intensity distribution in the lattice, /phi(2;) — Jo cos^ (nx/d), for d = 
7r/5, 2n, and IOtt, respectively. The two-dimensional intensity field in the respective solitons is shown in panels (d), (e), and (f), and 
the corresponding intensity profiles in two cross sections of the solitons, along y = and a; = 0, are displayed in panels (g), (h), and (i). 

[15] M. G. Vakhitov and A. A. Kolokolov, Izv. Vuz. Radiofiz. 16, 1020 (1973) [in Russian; English translation: Sov. J. Radiophys. 
Quantum Electr. 16, 783 (1973)]; see also L. Berge, Phys. Rep. 303, 260 (1998). 
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Fig. 2: The total intensity (power) (a) and peal< intensity (b) of families of two-dimensional soli- 
tons, vs. the absolute value of the propagation-constant shift, /i, for d = 2n, and different val- 
ues of the photonic-lattice strength Jo- The curves terminate at points where the solitons cease to exist. 
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Fig. 3: The minimum peak intensity, /p,min, and total power, Pmin, necessary for the existence of 
the two-dimensional solitons, vs. the strength, Iq, of the underlying quasi-one-dimensional lattice. 




Fig. 4: (a) Stable evolution of a high-intensity soliton, with Jo = 5 and Ip = 56.76, under random perturbations, is il- 
lustrated by the z-dependence of the soliton's peak intensity, (b) The same for a low-intensity soliton, with Ip — 1.1. 
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Fig. 5: The maximum value of the boost parameter, Q^, which generates solitons moving (actually, tilted 
in the spatial domain) along the x axis, vs. the soliton's peak intensity, Jp, for different value of the 
strength Jo of the underlying lattice. For Qx > Qmax, the application of the boost destroys the soliton. 




Fig. 6: A typical example of destructive collisions between moving solitons, for Jo = 5 and Ip = 1.78. The velocity vec- 
tors of the colliding solitons (shown by arrows in the first panel, in this figure and below) have components Qx.y = ±0.12. 
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FIG. 7: Merger of colliding solitons which were set in motion by the boost = ±0.01, 

with a finite aiming mismatch in the y direction. In this case, Jo = 25 and Ip — 83.97. 
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FIG. 8: Merger as a result of head-on collisions between slowly moving solitons that were set in motion (a) along 
the diagonal direction, by the application of the boost with Qx ~ Qy ~ ±0.01, and (b) in the vertical direction, 
by the boost Qy = ±0.01. In case (a), Jo = 25 and Ip = 83.97; in case (b), Jo = 15 and Ip = 27.92. 
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FIG. 9: Merger of broad (multi-lobed) solitons caused by the head-on collision. The solitons were 

set in motion by the vertical boost with Qy = ±0.01. In this case, Jo = 5 and Ip — 1.786. 



